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\^ I Abstract: Motivated by the fact that neutrinos are massive, we study the effect of 

neutrino Yukawa couphngs on neutrahno dark matter observables within the framework 
^^ ■ of a supersymmetric seesaw. We find that neutrino couphngs significantly affect the 

^T) • neutrahno rehc density in regions of parameter space where soft SUSY-breaking slepton 

■»H- ■ masses and/or trihnear couphngs are large. Depending on the size of the couplings, the 

^^ , neutrahno relic density spans over an order of magnitude in the A-funnel, focus point and 

f~^ I stop-coannihilation regions of mSUGRA. We also show that dark matter detection rates 

can be modified by up to several orders of magnitude. 
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1. Introduction 

A very attractive aspect of i?— parity conserving supersymmetry (SUSY) ||l|, Q is that it 
predicts the existence of a massive, electricaUy and color neutral, stable weakly interacting 
particle that is a natural Cold Dark Matter (CDM) candidate. In most cases this particle is 
the lightest neutralino, Z\, whose relic abundance can be reliably calculated as a function 
of model parameters ||3[ . Cosmological measurements determine the mass density of CDM 
with high accuracy. A combination of WMAP data with distance measurements from the 
baryon acoustic oscillations in galaxy power spectra gives Q 

9.h^ = 0.112010:00?^ (2a) , (1.1) 

where ri = pj pc with pc, the critical mass density of the Universe, and h is the scaled 
Hubble parameter. Such a precise measurement puts severe constraints on new physics 
scenarios. For example, in the well-studied minimal supergravity model, mSUGRA (or 
CMSSM) Q , the only surviving regions of parameter space are those in which neutralino 
annihilation is enhanced: the bulk region P, 1^, the stau [0, 0] or stop [10| coannihilation 



regions, the hyperbolic branch/focus point (HB/FP) region [11|, and the Aor h resonance 
annihilation (Higgs funnel) regions [0, |l2|, |l^ . The narrowness of these regions in mSUGRA 
motivated studies with one and two parameter extensions, in which non-universality of soft 



SUSY-breaking parameters reduce Vly h in accord with WMAP data over a large portion 



of parameter space ||14|. 



Other new physics that must be incorporated into the MSSM is that neutrinos are mas- 



sive, as indicated by the observation of neutrino ocillations [15|. An elegant explanation for 
neutrino mass is offered by the seesaw mechanism [^. Here three right-handed neutrinos 
(RHNs) with large Majorana masses are added to the SM. The decoupling of these RHNs 
naturally generates light masses for the left-handed neutrinos. We demand that the MSSM 
accommodate neutrino masses via the seesaw mechanism and explore whether favorable 
relic neutralino densities ensue by using the mSUGRA framework as a specific example. 
(Incidentally, this is potentially a more compelling scheme to expand the allowed mSUGRA 
parameter space than the introduction of non- universal soft SUSY-breaking parameters.) 
We demonstrate the invalidity of the lore that neutrino Yukawa couplings below the Grand 
Unification Scale (GUT) do not significantly affect Dark Matter (DM) observables.^ 

The rest of the paper is organized as follows. We briefly review the SUSY-seesaw 
mechanism in section ^. In section ^, we study the effect of neutrino Yukawa couplings 
on DM observables using benchmark points from mSUGRA augumented with RHN super- 
fields (mSUGRA+RHN). We show that the presence of large neutrino Yukawa couplings 
significanty affect the evolution of sparticle masses with concomitant changes in DM ob- 
servables. Finally, in section ^ we discuss the impact on the mSUGRA+RHN parameter 
space and on various non-universal models, and summarize our results. 

^Recently it was shown in the context of particular SUSY-GUT models that RHNs significantly impact 
low-energy phenomenology |l7|. However, the result is dominated by additional contributions to right- 
handed slepton mass evolution above the GUT scale, which is highly model dependent. In this work we 
take a model-independent approach and study the effect on DM predictions due to large neutrino Yukawa 
couplings below the GUT scale. 



2. SUSY-seesaw 

The superpotential for the MSSM with right-handed neutrinos (MSSM+RHN), in the 
notations and conventions of Ref . [^ , can be written as 

/ = Imssm + {UheabllHiN^ + ^(M;v).jiVfiV| , (2.1) 

where /mssm is the MSSM superpotential, L and Hu are, respectively, lepton doublet and 
up-higgs superfields, N[ is the superfield whose fermionic component is the left-handed 
antineutrino and scalar component is pjj, i,j = 1,2,3 are generation indices, a,b are 
SU{2)l doublet indices, i^ is the neutrino Yukawa coupling matrix and Mjq represents the 
(heavy) Majorana mass matrix for the right-handed neutrinos. After electroweak symmetry 
breaking, Eq. ( |2.lD leads to the the following mass matrix for the light neutrinos: 

M, = i.M^Xvl = ^vl , (2-2) 

where Vu is the vacuum expectation value of the up- type Higgs field h^. This corresponds 



to the so-called type I seesaw ||1^, where there are no other contributions (i.e., type II or 
type III) to the light neutrino masses. 

Additional soft SUSY breaking (SSB) terms are included so that the Lagrangian be- 
comes 



r,t /_2 



^ = ^MSSM - '^Rii^"^uii)ij'^Rj + 



.^r,t 1 



{aiu)ijeabLlH^i)\- + -{K)iji'Rii>Rj + h.c. 



(2.3) 



The parameters ^/(in? )ij, {sLu)ij and {hi^)ij are assumed to be of order the weak scale, 
even though right-handed neutrinos and their superpartners have much larger masses Mj^. 

Without further assumptions the neutrino Yukawa couplings are arbitrary. We can 
estimate the size of f,y by turning to SO(IO) GUTs. Here, all SM fermions and the right- 
handed neutrino of each generation are unified in a single spinorial representation, 16, of 
the SO{10) gauge group. The product of two 16's, that appear in the Yukawa term of a 
superpotential, can only couple to 10, 120 or 126. If higgs superfields reside in 10, as 
in the simplest models, then f^, = f^ at Mgut', if the higgses occupy 126, then (^ = Sfu- 
Contributions from 120 can only be subdominant, since they would lead to at least a pair 
of degenerate heavy up-quarks [^] . 

In the renormalization group equations (RGEs), the Yukawa matrix appears as fjfi^ and 
is dominated by the (3,3) element; if f^, = f^ then (fJf;/)33 ~ 0.25 at Mgut, while the other 
elements are smaller by two to five orders of magnitude. Also, the off-diagonal entries in 
the SSB matrices at the weak scale have to be very small to meet stringent fiavor-violation 
constraints and consequently they do not affect the mass spectrum significantly |19|| . 
Therefore, for the sake of clarity, we present our discussion in the third-generation dom- 
inant scheme, i.e., we assume f„,d,e,i/ ~ {fuAe,^)^^ = ft,b,TM, Mn ~ (^^Af)33 = ^^3, 
^u,d,e,u ~ {^u,d,e,u)^2 = —■^t,b,T,u ft,b,T,v and the SSB mass squared matrices are diagonal at 
all scales. However, the numerical analysis is performed in the full matrix form.^ 

^We assume a normal hierarchy for f^ and Mn for consistency with the other SM fermions. This does 
not restrict the light neutrino mass eigenstates to have a normal hierarchy. 



The right-handed neutrino superfields above the seesaw scale modify the evolution of 
gauge gi , 52 coupling at the 2-loop level, charged lepton and up-quark Yukawa couplings 
at the 1-loop level and introduce new RGEs for neutrino Yukawa couplings and Major ana 
masses [^0|. The RGEs for the parameter bjy are irrelevant for our analysis. Among the 
RGEs for the SSB parameters, the following get modified at the 1-loop level [ 
are relevant for our discussion: 
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where t = logQ, X„ 



m 



L-s 



m,r,_ + 171%, + A^, and Xt, X},, Xr, S and Cj, c" are given 



I^R 



H-n 



in Ref. |22|. The rest of the RGEs remain unchanged at the 1-loop level from those for the 



MSSM and may be found in Ref. |23|. Below the seesaw scale the right-handed neutrinos 
are integrated out and the MSSM is the effective theory with an appropriate change of 
RGEs; the only remaining trace is in the evolution of the dimension-5 neutrino operator 
k|2|. 

Usually, the f^Xt term dominates over gauge-gaugino terms and causes radiative elec- 
troweak symmetry breaking (REWSB) by driving m^ to negative values. In the case of 
the MSSM+RHN, mjj is driven to more negative values by the f^Xi, term. Similarly, the 
third generation slepton doublet mass m'j^ is driven to smaller values in the MSSM+RHN 
by the f^Xy term. Unless the trilinear A-terms have very large GUT-scale values, they 
are pushed to negative values by the gauge-gaugino terms. 

These features are illustrated in Fig. Q, which shows the running of third generation 
slepton doublet, up- and down-higgs mass parameters and A-terms from Mqut to My^^ak 
for mSUGRA (solid curves) and two cases of mSUGRA+RHN with different values of 
GUT-scale neutrino Yukawa coupling fi, (dashed and dotted curves). In the left frames 
we consider the case in which large Xy values result from large scalar masses. The trilinear 
A-terms evolve identically for the three models because the smallness of the A-terms at 
the GUT scale nullifies the upward push of fy. As expected, m,^ and m^ run to smaller 
values for mSUGRA+RHN. We see kinks in their evolution at the scale Q ~ Mj^.^, where 
the right-handed neutrino decouples and the f^Xy and fyAy terms do not contribute to 
the RGEs. Since right-handed neutrinos do not couple to the down-higgs directly, the RGE 
for m?^ is unaffected. In the right frames we consider the case in which large Xy values 
result from large trilinear A-terms. Now, the effect of fy in Eqs. ( 12.6^2.71 ) is non-negligible 
as is evident from the smaller absolute values for Aj- and At at the weak scale. As in the 
previous case, the large Xy pushes the slepton doublet mass to smaller values. However, 
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and 



''H l\l \"''H I' ''^^"= uiiiiu gciicicLuiuii Dicpuuii LLuuuicu iiidoc 1 1 1, j^ , aud triHuear couplings as a 
function of the renormalization scale Q in the mSUGRA model (solid), and two mSUGRA+RHN 
models - with fv^Mcur) — ftiMcur) (dashed hues) and with /^{Mgut) ~ '^ftiMcuT) (dotted 
lines). The model parameters are ttiq = 1507 GeV, mi/2 = 300 GeV, Aq — for the left frames 
and mo = 300 GeV, mi/2 = 500 GeV, Aq = -1200 GeV for the right frames. We take tan/3 = 10, 
^ > and mt = 171 GeV and also fix Mn^ = lO^^ GeV for mSUGRA+RHN models. Grey bands 
indicate the scale at which the third generation right-handed neutrino decouples. 



the picture is different for the up-higgs mass. Initiahy, the large Ai, contribution to X^^ 
pushes m'jj to more negative values. This reduces the magnitudes of Xi, and Xt, thus 
increasing the role of the gauge- gaugino terms. Consequently, weak-scale values of mj^ 
are almost unchanged from those in mSUGRA. 

For moderate to large tan/3 values, as favored by the Higgs boson mass constraint 
from LEP2, the tree-level^ minimization condition for EWSB in the MSSM can be ap- 
proximated as /U^ ~ ~'^H ~ 0.5M^, and the the tree-level CP-odd Higgs boson mass is 
m\ ~ m?^ — rnjj — M^, which implies that when [m^ | is large, both |/i| and niA are 
large and correspondingly the H and H Higgs bosons are heavier. Also note that larger 
weak-scale A-terms lead to more mixed states and a lighter /i sfermion state. 



^We use tree-level relations only for illustration. In practice, fi^ and sparticle spectrum are determined 
using the full 1-loop expressions. 



3. Seesaw and Relic Density 



In the MSSM, the neutrahno is a superposition of the bino, wino, up- and down-higgsino 
states. Over the most of the mSUGRA parameter space the hghtest neutrahno Z\ is mainly 
bino with a smaU annihilation cross section so that its relic density is significantly above 
the WMAP range. Nevertheless, there exist regions of parameter space where various 
mechanisms enhance the annihilation and WMAP bounds can be satisfied. To exam- 
ine effects of a SUSY-seesaw on DM observables in these regions, we use the mSUGRA 
model as an example. The model is completely specified by the well-known parameter set 
"T-O) "ii/2i ^Oi tan/3, sgn(//). The SSB terms are taken to be universal at the GUT scale, 
^GUT ~ 2 X 10^^ GeV, with uiq the common scalar mass, m-i /2 the common gaugino mass, 
and ^0 the common trilinear coupling. The universal SSB boundary conditions allow us 
to isolate the neutrino Yukawa coupling effects. Since electroweak symmetry is broken 
radiatively, the magnitude (but not the sign) of the superpotential Higgs mass term /i is 
determined, and we can trade the GUT-scale bilinear soft term B for the weak-scale ratio 
of Higgs vacuum expectation values, tan/?. Once weak-scale values of the SSB parameters 
are computed via renormalization group evolution, they serve as inputs for computation of 
sparticle masses and mixings. Then, the neutrahno relic density Q^ h? and DM detection 
rates can be calculated. 

In each DM-allowed region of mSUGRA we select one benchmark point "^ from Table 
and vary MjVg . The values of the first and second generation right-handed neutrino masses 
are not important, since the corresponding Yukawa couplings are assumed to be very 
small. For definiteness, we fix their masses to 10^^ GeV and 10^^ GeV, respectively. For 
the GUT-scale neutrino Yukawa coupling we consider the two cases introduced earlier: 
MMgut) = MMgut) and MMqut) = SftiMcur)- We take ^ > as suggested by 



experimental measurements of the muon anomalous magnetic moment |25], and top mass 
rrit = 171 GeV to conform with Tevatron data [p^. 

To properly account for RHNs, 
we upgraded ISAJET ^] to full 
matrix form and added RGEs for 
evolution of the RHN Majorana mass 
matrix Mjv above the scale of de- 
coupling, and a dimension-5 opera- 
tor K below it |28|. To evaluate the 
neutrahno relic density^ and direct 
DM detection rates we employ the 
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Table 1: Input parameters for benchmark points and 
corresponding DM-allowed regions of mSUGRA. For all 
points /i > and nit = 171 GeV. 



IsaRED g and IsaRES |29|] subrou 

tines of the IsaTools package. To 

compute the indirect DM detection rates we use the DarkSUSY 



package. We evalu- 



The small difference in parameter values 



''These points are closely related to SPS benchmark points 
is attributable to the different code and top mass we are using. 

^We adopt the conceptually simplest scenario in which neutralinos are thermally produced in the stan- 
dard ACDM cosmology and make up the DM. 
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Figure 2: Neutralino relic density as a function of the GUT-scale value of M^^ for two values of 
the GUT-scale neutrino Yukawa coupling: f,y{McuT) = ftiMcur) and f„{McuT) = 5ft{McuT)- 
The labels on the curves refer to the points in Table and the green bands mark the WMAP range 
of Eq. (1.1). The curves do not extend below the value of A/jy., for which the tau- neutrino becomes 
heavier than 0.3 eV. 



ate the integrated continuum 7-ray flux above a E^ = 1 GeV threshold. For positrons 
and antiprotons, we compute the differential flux at a kinetic energy of 20 GeV, for 
which optimal statistics and signal-to-background ratio are expected at spaceborne an- 



tiparticle detectors |31]. For antideuterons, we compute the average differential flux in the 
0.1 < Tfj < 0.25 GeV range, where Tg is the antideuteron kinetic energy per nucleon P2|| . 



We performed cross-checks against ISAJET 7.77 with IsaTools as well as SPheno 2.2.3 |33| 



interfaced with micrOMEGAs 2.0.7 |34] and found results in good agreement. 

The effect of neutrino Yukawa couplings on the neutralino relic density is illustrated 
in Fig. |2[ which shows O^ h"^ versus the GUT-scale value M^^ for the points in Table |l|. 
The curves do not extend below the value of M^^ for which the tau-neutrino^ becomes 
heavier than 0.3 eV, in accordance with the cosmological bound on neutrino masses, 
Yl^^f ^ 1 sV IQ]. The curves become flat for MjVg ^ 2 x 10^^ GeV because the right- 
handed neutrino decouples above the GUT scale and has no effect on the spectrum. We 
see that, contrary to common belief, neutrino Yukawa couplings significantly impact the Zi 
relic density. In what follows, we provide detailed explanations of this effect and discuss 
consequences for DM detection rates. 

3.1 Bulk region 

If the neutralino is bino-like, as in many SUSY models, and sparticles are light (~ 100 GeV) 
then Zi can efficiently annihilate into a pair of SM fermions via a sfermion exchange in 
the t— channel. The dominant process is ZiZi — > U, since right-handed sleptons have the 



® Since flavor effects are irrelevant for our study, we ignore mixing in both quark and lepton sectors, thus 
treating neutrino favor eigenstates as mass eigenstates. 
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Figure 3: Some relevant sparticle masses (upper frames) and enhancement factors, i.e., ratio of 
the rate in mSUGRA+RHN to the corresponding rate in mSUGRA, for direct and indirect DM 
detection (lower frames) as a function of Mns at the GUT scale for the benchmark point in the 
bulk region. The legend is described in the text. The neutrino Yukawa coupling is f^iMcuT) — 
ftiMcur) in the left frames and fv{McuT) = ^ftiMcur) in the right frames. 



largest hypercharge value (Y{E'^) = 1). In mSUGRA this occurs in the bulk region of 
parameter space which has small values of mo, m-^ /2, ^o- This region was favored in early 
work on mSUGRA, but is largely excluded by the nonobservance of chargino and slepton 
pair production at LEP2 ||6|, 0]. 

In the upper frames of Fig. ^, we plot relevant sparticle masses and the /x parameter 
versus the GUT-scale value Mjvg for point A of Table |l|. We see that the tau-sneutrino gets 
lighter as Mn^ decreases - a smaller M7V3 means a lower seesaw scale and a correspondingly 
greater downward push from the f^Xy term in Eq. (p.4|). However, due to the smallness 
of m-o and ^0, the overall effect is small - up to ~ 3% for fuiMcuT) = fti^GUT) and 
up to ~ 15% for fuiMcuT) = '^ftiMcuT)- Because of the universality of GUT-scale SSB 
boundary conditions, the lightest stau is mostly right-handed and hence its mass is affected 
only slightly. On the other hand, f2 is mostly left-handed and experiences changes similarly 
to Dt- The f^A^, term in Eq. ( |2.7D makes At less negative leading to a slightly heavier ii: 



by up to ~ 1% in the left frame and by up to ~ 3.5% in the right frame. The rest of the 
masses are unaffected by the right-handed neutrino. As a result the neutralino relic density 
is barely affected; see curves A in Fig. ||. 

The lower frames of Fig. ^ show the effect of changing the right-handed neutrino 
mass on neutralino DM direct (DD) and indirect detection (IDD) rates. It is well-known 
that the rates are sensitive to the (unknown) DM halo distribution as well as values of 
hadronic parameters ^ ^. We define an enhancement factor as the ratio of the rate 
in mSUGRA+RHN to the corresponding rate in mSUGRA, so that it is approximately 
independent of those uncertainties. We see that as Mjy^ decreases, the rates for positrons 
(orange dashed curve), muons from the Sun (blue dash-dotted curve) and 7-rays from the 
Galactic center (green long-dashed curve) increase. This is because as M7V3 decreases, -f2 
gets lighter, and neutralino annihilation in r lepton pairs is enhanced with a concomitant 
enhancement of the above IDD rates. The difference in enhancement factors corresponds 
to the role r leptons play in these rates: r's are the dominant source for positrons and 
muons (since vector boson production is kinematically forbidden for the parameter point 
considered), but contribute little to gamma-ray production Q. For the fv{MGUT) = 
ft{MGUT) case, rates for positrons, muons and 7-rays increase by up to 30%, 20% and 
2% respectively. In the case of fv{MGUT) = ^fti^GUT), the f2 mass is pushed closer 
to the Zi mass causing larger enhancements of the IDD rates. Rates increase by up to 
160% for positrons, 80% for muons and 20% for 7-rays as shown in the lower right frame 
of Fig. y. Direct detection rates, that are conventionally represented as the cross-section 
of spin- independent elastic neutralino-proton scattering (solid magenta curve), as well as 
the antiproton (dashed black curve) and antideuteron (solid turquoise curve) rates are 
unaffected by the ZiZi — > rf enhancement; their tiny suppression (less than 1%) is due to 
the small change in sbottom mixing caused by radiative corrections. 

3.2 Stau-coannihilation region 

If the mass of a bino-like neutralino is close to the mass of a more strongly interacting 
sparticle, then fi^ h? is lowered via coannihilation. Usually the lightest stau fi acts as 
such a sparticle and rapid reactions Zifi -^ X and fifi — > X (where X denotes any 
allowed final state of SM or Higgs particles), in the early universe make the Zi relic density 
consistent with WMAP data. In mSUGRA this region appears at small niQ and small to 
medium ?n-i/2 values |§, |9[. 

The situation is similar to that in the bulk region - Dj- and f2 get lighter as MJV3 
decreases, while the rest of the sparticles are unaffected. However, the diminution is smaller 
(< 1%) because Aq = and mi/2 > "^o reduces the relative role of the fi, terms in the 
RGE evolution. Again, due to the right-handedness of fi, its mass and concomitantly 
n^ h? remain unchanged - see curves B in Fig. |2|. 

The smaller change in the f2 mass leads to a smaller increase in the ZiZi -^ tt 
rate, which in turn yields a smaller enhancement for muon and positron IDD rates (< 2% 
for fy = ft at the GUT scale and < 12% for f^ = 3ft) as compared to the bulk re- 
gion. Also, in this region, neutralino pair annihilation into vector bosons is allowed since 



m^ ~ 0.5?7ii/2 > Mz- This process dominates over Zi pair annihilation into r leptons as 
the source of muons and positrons and is unaffected by the RHN mass; this further reduces 
the sensitivity of these IDD rates to changes in the f2 mass. The effect of varying M7V3 on 
the DD and the other IDD rates is the same as in the bulk region and is of the same order. 



3.3 Stop-coannihilation region 

Another possibility for coannihilation is with the scalar top of mass mi ~ tti^ . In 
mSUGRA this region is located at small mo and raii2 and is characterized by a large 
value of ^0 0- 

For l^ol ^ "^0) ?^i/2; the A-terms dominate the RGE evolution. The large negative 
f^Ay term pushes \At\ and \Ar\ to smaller values at the weak scale, which reduces the L-R 
mixing in stops and staus. The smaller \At\ contribution to Xt also reduces the downward 
push by the ffXt term in m? evolution. A combination of these effects increases mr away 
from the neutralino mass, thus shutting off the coannihilation mechanism. Also, the large 
Xy term in Eq. (2.4) drives the masses of the left-handed sleptons Vr and f2 to smaller 



values; see Fig. |^. For fy{MGUT) = fti^GUT)-, the growing stop mass causes the neutralino 
relic density to rapidly increase above the WMAP-preferred range to O^ h?' ~ 0.3 for 
^Na < 2 X 10^^ GeV; see curve C in the left frame of Fig. |. 

The dependence on Mn^ is more significant for the fu{MGUT) = "iftiMcuT) case, 
shown in the upper right frame of Fig. ^. Such a large neutrino Yukawa coupling increases 
the importance of the fy terms in the RGEs. The stop mass is pushed to even higher values, 
while m? falls faster and causes stau masses to decrease by ~ 23%; staus also undergo an 
identity flip at M7V3 — 10^^ GeV - fi changes from a right-handed to dominantly a left- 
handed state. But most importantly, the tau-sneutrino gets lighter (with mass ~ ttilz) and 
becomes the NLSP with almost the same mass as the neutralino for M7V3 ^ 4.5 x 10^^ GeV. 
This opens up the Zi — Vj- coannihilation channels that rapidly decrease Q.^ h? down to 
and then below the WMAP range, as shown by curve C in the left frame of Fig. ^. Below 
M7V3 — 3 X 10^^ GeV, the u^- gets lighter than the Zi and the parameter space closes. 

As M7V3 is dialed to smaller values, tau-slepton masses get smaller, leading to enhanced 
Zi pair annihilation into r-leptons in the t— channel. Since r's are one of the primary 
sources for muons and positrons, the corresponding IDD rates increase by up to 100% and 
150% respectively in the fuiMcuT) = fti^GUT) case. Similarly, 7-rays originate from r's 
and therefore their IDD rate changes by up to 130%. In the fuiMcuT) = 'ifti^GUT) case, 
significantly lighter staus cause the IDD rates to increase by up to 300% — 700%. 

Antiprotons and antideuterons are mainly produced in the hadronization of heavy 
quark flavors (bottom and top) present in the products of Zi pair annihilation. In general, 
increasing the ti mass suppresses annihilation rates into ti. However, since neutralinos in 
the halo are essentially at rest, the process ZiZi -^ ti is kinematically forbidden for the 
point considered. Sbottoms, on the other hand, are very heavy and do not significantly 
contribute to ZiZi — >■ 66. As a result, IDD rates for antiprotons and antideuterons change 
by less than 1% as M7V3 is varied. 
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Figure 4: Similar to Fig. a but for the benchmark point in the stop-coannihilation region. 



3.4 Higgs funnel 



If 2m^^ 



rriA, neutralino annihilation through the A (and H) Higgs boson in the s— channel 
are resonantly enhanced. Since the ^- width can be very large (F^ ~ 10 — 50 GeV), an exact 
equality in the mass relation is not necessary to achieve the near-resonance enhancement. 
This WMAP-compatible region is known as the ^-funnel fl^ and occurs in mSUGRA at 
medium uiq and large tan/?. 

Increasing rriQ increases the downward push by the fyX^ term in Eq. (2.4), leading 
to a lighter Vj- and f2. Since A^ = 0, the L-R mixing is small, fi is right-handed and 
therefore almost unaffected. Similarly, the f'^Xy term in Eq. (|2.8| ) pushes rn^j to more 
negative values resulting in larger |//| and mA values. This is borne out in the upper 
frames of Fig. ^. Larger |/i| translates into heavier (< 7%) higgsino-like Z3, Z4, W2 states. 
Increasing the mass of the CP-odd Higgs boson A pushes it away from the resonance 
resulting in a reduction of the Zi annihilation rate and thus larger relic density. From 
curves D in Fig. ||, we see that fi^^/i^ is above the WMAP range for M^^ < 10^^ GeV (left 
frame); a larger neutrino Yukawa coupling causes rriA to grow faster and il^ /i^ exceeds 
the WMAP range for Mn^ < 10^*5 GeV (right frame). 



10 



(0 
(0 

to 

s 



600 
550 
500 
450 
400 
350 
300 
250 
200 
150 



c 
q> 

E 

u 

c 

c 



point D: nig = SOOGeV, my2 = 450GeV, ten(3 = 51, Ag =0 

fv(Mn„T) = f,(MnirT) 



1 1 1 1 M 1 1 


' L' L" L" L' L' L' - ^' L' U U U U U U. « M. « n F^ Fh Fh Fh n Fh n 


^2 

A 

^1 


E ^ 








I ^\ 




^ v 




- "1 





1 


- 






^^. — 


0.8 


- 


-"^ ^^^^ 


:^" 




0.6 


— 


^^^^^ 














DDs/ 


0.4 
0.2 


;i 




- • ^^ 


muons f/SuA7 

gammas 

positrons 

antiprotons 

antideutrons 



10 



13 



10 



14 



10 



15 



16 



17 



10 10 





fv(M^UT) = 3f,(M^UT) 




L ^i 


" " !_' ^' L,' »■■ f . 1^ 1-. •■. 1-. r. r. r. 1- 


^2 

A 




E V, 


_,_----^^*<--.....^ 




= ^^1 




— 7 




z "1 





/f 





yyi 




X' 




^y 




^^ y / 




^---^ / / 




L ^y ■:- 


DDsi 

muons f/Sun 

gammas 

positrons 

antiprotons 

antideutrons 



10 



15 



10 



16 



io'' 



Figure 5: Similar to Fig. ||, but for the benchmark point in the A-funneL 

The same neutrahno annihilation mechanism ZiZi — > A — > bb/rf that reduces the 
reUc density plays a primary role in DM halo annihilation. Consequently, its reduction has 
a significant effect on IDD rates. Increasing the mass of the A decreases the production 
of 5-quarks in the Zi halo, thus reducing antiproton and antideuteron fluxes; the flux of 
7-rays produced via the 6 — > vr^ ^ 7 chain also experiences similar suppression. From the 
lower frames of Fig. ^ we see that the effect can reach ~ 40% if fv{MGUT) = ftiMcuT) 
and ~ 75% if U{Mgut) = SftiMouT). 

The exchange of h and H bosons are important for neutralino-proton elastic scattering. 
Therefore, increasing mn leads to a suppression of DD rates by up to about 15% in the left 
frame and up to ^ 50% in the right frame. The spin-dependent Zi—p elastic cross section 
suffers a similar suppression thus lowering the neutralino capture rate in the Sun. This 
effect is coupled with a reduction of the neutralino annihilation cross section. Together, 
the solar muon flux is reduced by ~ 20% for fu{MGUT) = fti^GUT) and by ~ 75% for 
UMGUT)=^ft{MGUT)- 

3.5 Hyperbolic branch/Focus point region 

If /i is sufficiently small the lightest neutralino develops a sizable higgsino component, which 
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Figure 6: Similar to Fig. y, but for tlic benchmark point in tlic HB/FP region. 

enhances Zi annihilation to WW, ZZ and Zh in the early universe and brings O^ h? in 
accord with Eq. (p..lj). In niSUGRA this is realized in the hyperbolic branch or focus point 
region located at very large rriQ along the edge of the no-REWSB region where ^j? becomes 
negative |11]. 

The presence of a large neutrino Yukawa coupling in Eq. ( |2.8| ) results in a larger 
absolute value for m?^ at the weak scale and a corresponding increase in /i. This means 
that the neutralino becomes increasingly bino-like and its relic density rapidly grows beyond 
the WMAP range, as shown by curves E in Fig. ||. The lightest neutralino mass increases 
by about 15%, while the lightest chargino mass grows by 50% and changes composition 
from higgsino to wino; see the upper frames of Fig. ^. The masses of heavier charginos and 
neutralinos also increase by 50% to 100%, as \ii\ increases. As |?7z^ | increases, the A-mass 
also increases by up to 2.5% (up to 10% for the larger f^ case). As in Section p.3|, decreasing 



m? causes the mass of the Vr to drop by ~ 3% (up to ~ 15% for the larger fy case) and the 
fi to change its composition to left-handed state at Mj^.^ < 2 x 10^^ GeV {M^.^^ < 10^^ GeV 
for the larger fi, case) with an accompanying reduction in mass. However, due to the large 
mo value, all sfermions are very heavy and are not relevant for DM detection. 
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The decrease of the higgsino composition of Z\ has a significant effect on DM detection 
- as M7V3 decreases, all the rates drop by several orders of magnitude. The lower higgsino 
content reduces annihilation into vector bosons and suppresses the overall Zi annihilation 
rate. As a result 7-ray and antimatter fluxes from neutralino annihilations in the galactic 
halo fall by a factor of 10^ (~ 10^ in the right frame). DD rates decrease by up to about two 
(three) orders of magnitude because the coupling of Zi to higgs bosons is diminished. The 
muon flux is also reduced due to a decreasing neutralino capture rate and weakening Zi 
pair annihilation: the overall reduction reaches two orders of magnitude for fuiMcuT) = 
ftiMcuT) and more than six orders of magnitude for fuiMcuT) = ^ft{MGUT)- 

4. Discussion and Conclusions 

So far our discussion was confined to particular values of tuq, ""1-1/2) ^O) but neutrino 
Yukawa couplings generally affect the WMAP-allowed regions in mSUGRA parameter 
space. To illustrate this we performed a random scan in (mo, m-1/2, M^^) and plot- 
ted results in the conventional (jit-Oi "^1/2) plane after marginalizing over Mi^^. In the 
green regions of Fig. 0, the neutralino relic density lies within the WMAP 2(T-range given 



by Eq. (1.1). We do not show points that are excluded by the LEP2 chargino bound 
m^ > 103.5 GeV |^^. For comparison, we superimposed the WMAP-allowed regions in 
mSUGRA as blue crosses. The stau-coannihilation region is virtually unaffected by the 



neutrino Yukawa coupling - as explained in Section 3.2. On the other hand, the fi param- 
eter in the HB/FP region gets larger for smaller MiVg, thus postponing the breakdown of 
the REWSB mechanism to larger itiq values. This makes the HB/FP region very sensitive 
to the neutrino Yukawa coupling - depending on the value of M^^ it can move to the right 
by up to ~ 500 GeV, and even more if fi, = 3 ft at the GUT scale. This shift of the HB/FP 
region to larger mg values makes sfermions heavier and can have a sizable effect on collider 
signatures. 

If the mechanism that lowers the relic density also affects DM detection rates, then rates 
in the WMAP allowed regions may be qualitatively unaltered from mSUGRA expectations, 
as is the case in the HB/FP region where both the relic density and detection rates are 
tied to the higgsino component of Zi. This is because neutrino Yukawa coupling effects 
can be compensated by moderate shifts in the space of fundamental parameters. However, 
when the mechanisms dictating the relic density and detection rates are different, and such 
a compensation is not possible, one can expect sizeble variations in DM detection rates. 
In the bulk and stau-coannihilation regions, the relic density remains unaffected by dialing 



M7V3 , but the rates can change appreciably - see Sections |3^ and |3.2| . Since the specific 
realization of the seesaw mechanism is unknown, it is an additional uncertainty in the DM 
detection rates. 

Although we presented our analysis of the effects of neutrino Yukawa couplings in a 
mSUGRA+RHN framework, the effects appear in any SUSY-seesaw model in which RHNs 
with large Yukawa couplings decouple below the GUT scale. Moreover, the effect can be 



even larger in a non-mSUGRA framework due to different sparticle properties [14|. For 



instance, the right-handedness of the fi that limited the effect of neutrinos in mSUGRA 
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Figure 7: Regions compatible with the dark matter relic density in the (mo, f^i/2) plane for 
the mSUGRA+RHN model with U{Mgut) = MMgut), Aq = 0, fj. > and rut = 171 GeV. 
The green points have a neutralino relic density within the WMAP 2cr-range (Eq. ^]|). The right- 
handed neutrino mass M^^ is randomly varied between 10^'^ GeV and Mgut — 2 x 10^^ GeV, 
with the requirement that tau-neutrino be lighter than 0.3 eV. The blue crosses mark the WMAP- 
allowed regions in the mSUGRA model. Note the significant expansion of the parameter space 
when neutrino masses are accounted for. 



(especially in the stau-coannihilation region) does not hold in scenarios in which one or both 
higgs mass parameters are non- universal at the GUT scale, as in the models of Ref. |p7|| . 
Here, fi is dominantly f^ and thus more susceptible to neutrino Yukawa couplings. Also, 
the ^-funnel appears in mSUGRA only at tan/3 where bottom and tau Yukawa couplings 
are large, thus dampening the effect of f^ on RGE evolution. In the models of Ref. p7|| , 
the ^-funnel can also occur at small and medium tan/3 values and we expect even larger 
changes of DM detection rates with variations of the RHN mass. 

We emphasize that although we used 50(10) models to estimate the size of the neu- 
trino Yukawa coupling at the GUT scale, our results are not limited to models that employ 
this gauge group. In fact, in S'[/(5)-based SUSY GUTs, where the right-handed neu- 
trinos are singlets, iy is not correlated with the other Yukawas and can be even larger 
than we considered. Requiring perturbativity, neutrino Yukawa couplings are limited to 
f. < 15fu 



It is also worth keeping in mind that many mechanisms other than the type I seesaw 
have been suggested for the generation of neutrino masses. For example, in a double 
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seesaw |16], one postulates an 50(10) singlet with mass Mg in addition to a right-handed 



neutrino A^, and finds that the light neutrino mass takes the form my ~ Msv^f^/Mf^ 



Since Ms <C M^^ , M^^ would be considerably smaller than in a type I seesaw: if Ms ~ 
1 TeV, the RHN could be as light as 10*^ GeV. In this case, changes in DM observables 
could be significantly larger. 

In summary, motivated by the concrete evidence that neutrinos are massive, we ex- 
amined the effect of neutrino Yukawa couplings on the neutralino relic density and dark 
matter detection rates in SUSY seesaw models. We found that, contrary to common be- 
lief, neutrino Yukawa couplings can significantly affect the relic density. Effects are most 
prominent in regions of parameter space where SUSY-breaking slepton masses and/or tri- 
linear couplings are large. Such conditions are satisfied, for example, in the j4-funnel, focus 
point and stop-coannihilation regions of mSUGRA. Here the neutralino relic density can 
change by over an order of magnitude due to the effect of neutrino Yukawa couplings; see 
Fig. ^. We also showed that DM detection rates can be changed by up to several orders 
of magnitude in the focus point region and by factors of a few to ten in the other regions; 
see Figs. pHBI. 
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